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ABSTRACT 

We show that detecting real roots for honestly n-variate 
(n + 2)-nomials (with integer exponents and coefficients) 
can be done in time polynomial in the sparse encoding for 
any fixed n. The best previous complexity bounds were 
exponential in the sparse encoding, even for n fixed. We 
then give a characterization of those functions k(n) such 
that the complexity of detecting real roots for n-variate 
(n + fc(n))-nomials transitions from P to NP-hardness as 
n — > oo. Our proofs follow in large part from a new com- 
plexity threshold for deciding the vanishing of ,4-discriminants 
of n-variate (n+fc(n))-nomials. Diophantine approximation, 
through linear forms in logarithms, also arises as a key tool. 
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1. INTRODUCTION AND MAIN RESULTS 

Consider real feasibility: the problem of deciding the exis- 
tence of real roots for systems of polynomial equations. In 
addition to having numerous practical applications (see, e.g., 
[BG-V03] !. real feasibility is an important motivation behind 
effectivity estimates for the Real Nullstellensatz (e.g., |Ste741 
SchOO ), the quantitative study of sums of squares [Ble041 
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IRS09I BHPR09 , and their connection to semi-definite pro- 
gramming and optimization Par03, Las07 . In particular, 
real solving of sparse polynomial systems arises in concrete 
applications such as satellite orbit mechanics [AM09] , and 
real solving clearly involves real feasibility as an initial step. 
We are thus inspired to derive new algorithms and complex- 
ity lower bounds for real feasibility, in the refined setting of 
sparse polynomials. 

To state our results, let us first clarify some basic notation 
concerning sparse polynomials and some well-known com- 
plexity classes. Recall that R* is the multiplicative group of 
nonzero elements in any ring R. 

Definition 1.1. When a } e 

(^l,j , • • • , O n ,j )? J = X^ ,J ' ■ ■ X n 



the notations aj — 
and x — (xi, . . . , x n ) 
where Cj G K* 



will be understood. If f(x) := Y7j=i CiXa 
for all j, and the aj are pair-wise distinct, then we call f 
a (real) n-variate m-nomial, and we define Supp(/) := 
{01, . . . , a m } to be the support of f. We also let T n , m de- 
note the set of all n-variate [m\ -nomials within 7J\x\ , . . . , x n ] . 
Finally, for any m > n + 1, we let J-n, m ^ J~n,m denote the 
subset consisting of those f with Supp(/) not contained in 
any (n — l)-flat. We also call any f G T^.m an honest 
n-variate m-nomial (or honestly n-variate ). o 

For example, l+lx'ix2x\x\—^x^ s 'x^x% %i xf 7 is a 4-variate 
trinomial with support contained in a line segment, but it 
has a real root i£l 4 iff the honestly univariate trinomial 
1 + 7yi — 43yi 9 has a real root y\ £ R. More generally 
(via Lemma 12.51 of Section [2.21 below) . it will be natural to 
restrict to J~n, n +k (with k > 1) to study the role of sparsity 
in algorithmic complexity over the real numbers. 

We will work with some well-known complexity classes 
from the classical Turing model, briefly reviewed in the Ap- 
pendix. (A more complete introduction can be found in 
Pap95].) In particular, our underlying notion of input size 
is clarified in Definition 12.11 of Section 12.11 below, and il- 
lustrated in Example 11.41 immediately following our first 
main theorem. So for now, let us just recall the basic inclu- 
sions NC 1 CPC NP C PSPACE. While it is known that 
NC 1 ^ PSPACE the properness of each of the remaining 
inclusions above is a famous open problem. 

1.1 Sparse Real Feasibility and ^-Discriminant 
Complexity 

Definition 1.2. Let R+ denote the positive real num- 
bers and let FEASr (resp. FEAS+j denote the problem 



of deciding whether an arbitrary system of equations from 
U n Z[xi, . . . , x n ] has a real root (resp. a root with all co- 
ordinates positive). Also, for any collection T of tuples cho- 
sen from \J k nef} (Z[x 1 , . . . ,x n ]) k , we let FEAS R (.F) (resp. 
FEAS+(jF) ) denote the natural restriction o/FEASr (resp. 
FEAS+ ) to inputs in T. o 

It has b een known since the 1980s that FEASr G PSPACE 
[Can88] , and an NP-hardness lower bound was certainly 
known earlier. However, no sharper bounds in terms of spar- 
sity were known earlier in the Turing model until our first 
main theorem. 

Theorem 1.3. Let Z + (f) denote the zero set of f in R". 
Then: 

0. FEAS+(U neN ^,„ +1 ) and FEAS R (U„ 6M ^,»+i) 
are in NC 1 . In particular, when f G Tn, n +\, Z+(f) 
is either empty or diffeotopi^ to R^ _1 , with each case 
actually occuring. 

1. For any fixed n, FEAS + (f* i „ +2 ) and FEASe(^,„+2) 
are in P. 



2. For any fixed e>Q, both FEAS+ (U„cn ^~n,n+n^, 
FEAS K (U neN ^"n,n+n=) are NP-hard. 



and 



Slightly sharper algorithmic complexity bounds hold when 
we instead work in the BSS model over R (thus counting 
arithmetic operations instead of bit operations), and this is 
detailed in [PRT09] . 

Example 1.4. A very special case of Assertion (1) of 
Theorem ] 1.3\ imvlies that one can decide — for any nonzero 
ci, . . . , C5 £ Z and D G N — whether 

ci + C2X? 9 + c z x\ 3 xf + c 4 xf D + c 5 xl 4 x?x 3 
has a root in R 3 , using a number of bit operations polyno- 
mial in 

tog(£>) + ]og[(|ci| + l)-.-(|cB| + l)]. 
The best previous results (e.g., via the critical points method, 
infinitesimals, and rational univariate reduction, as detailed 
in \BPR06f ) would yield a bound polynomial in 
D + log [(|ci| + 1) • • ■ (|c 5 | + 1)] instead, o 

We thus see that for sparse polynomials, large degree can 
be far less of a complexity bottleneck over R than over C. 
Theorem 1 1.3 1 is proved in Section [3 . 2 1 below. The underlying 
techniques include .4-discriminants (a.k.a. sparse discrimi- 
nants) (cf. Section [2.3[) . Viro's Theorem from toric geome- 
try (see the Appendix, or GKZ94, Thm. 5.6]), and effective 
estimates on linear forms in logarithms [Bak771 [Nes03j . 

In particular, for any collection Ta of n-variate m-nomials 
with support A, there is a polynomial Aa in the coefficients 
(d) called the ^-discriminant. Its real zero set partitions 
Ta into chambers (connected components of the comple- 
ment) on which the zero set of an / G Ta has constant 
topological type. A toric deformation argument employing 
Viro's Theorem enables us to decide whether a given cham- 
ber consists of / having empty or non-empty Z+(f). For 
any A C Z n of cardinality n + 2 (in sufficiently general posi- 
tion), there is then a compact formula for the „4-discriminant 
that enables us to pick out which chamber contains a given 
/: one simply computes the sign of a linear combination of 



logarithms. Our resulting algorithms are thus quite imple- 
mentable, requiring only fast approximation of logarithms 
and some basic triangulation combinatorics for Supp(/). 

Example 1.5. Consider A:= {(0,0,0), (999,0,0), (73,0, 19), 
(0, 2009, 0), (74, 293, 1)}, which gives us the family of trivari- 
ate pentanomials 

T A := {ci + c 2 x'f 9 + c 3 x?xf + Ci xj 009 + c 5 x\ i xf 3 x 3 \ c t Gl*} . 
Suppose further that f G Ta is an element satisfying ci , C2 , C3 , 
C4 > and C5 < 0. It then turns out via Lemma \2.13\ (cf. 
Section \2.J\ below) that Z+(f) has a degeneracy iff the A- 
discriminant, Aa{c) := 

-27886408 27886408 2677997 2677997 2006991 2006991 5561433 45561433 cf 132829 
vanishes. In fact, via the techniques underlying Theorem 
\1.S\ Z+(f) is either empty, a point, or isotopic to a 2-sphere, 
according as Aa(c) is positive, zero, or negative. Note in 
particular that determining the sign of Aa(c) is equivalent 
to determining the sign of 

38132829 log(3S132829) + 27886408 log(ci) + 2677997 log(c 2 ) + 2006991 log(c 3 ) + 5561433 log(c 4 ) 
-278864081og(2788S408) - 26779971og(2677997) - 2006991 Iog(2006991) - 5561433 log(5561433) - 381328291og(c 5 ). O 

While we review ,4-discriminants in Section [2.31 below, it 
is important to observe now how the computational com- 
plexity of ,4-discriminants closely parallels that of FEASpj 
compare Theorem [L3] above with Theorem [L7] below. 

Definition 1.6. Let ADISC= (resp. ADISC >y ) denote 
the problem of deciding whether A^(/) vanishes (resp. de- 
termining the sign of AA{f)) for an input polynomial f 
with integer coefficients, where A = Supp(/). Finally, let 
ADISC = (T) (resp. ADISC>(J' r ) y ) be the natural restric- 
tion o/ADISC= (resp. ADISC> ) to inputs in some family 
T. o 

Theorem 1.7. 

1. ADISC= I IJ Tn, n +2 ) £ P and, for any fixed n, 

VngN ' / 

ADISC>(^, n+2 )GP. 

2. For any fixed e > 0, both ADISC= I U JF* n+n e ) and 

ADISC> ( 1J TZ >n+n e ) are TSP-hard. 



1 Sec Definition [2~TU1 of Section [231 below. 



Theorem 1 1.71 is proved in Section [3.11 after the development 
of some necessary theory in Section [5] below. 

1.2 Related Work 

Earlier work on algorithmic fewnomial theory has 
mainly gone in directions other than polynomial-time al- 
gorithms. For example, Gabrielov and Vorobjov have given 
singly exponential time algorithms for weak stratifications of 
scmi-Pfaffian sets [GV04| — data from which one can com- 
pute homology groups of real zero sets of a class of functions 
more general than sparse polynomials. Our approach thus 
highlights a subproblem where faster and simpler algorithms 
are possible. 

Focussing on feasibility, other than the elementary results 
FEASi(fi (1 ) G NC° and FEAS,(fi, 2 ) G NC°, there ap- 
pear to have been no earlier complexity upper bounds of 
the form FEASr (^i, m )eP, or even F EASk (fi, ro )gNP, 
for m > 3. (With the exception of RY05 , algorithmic 
work on univariate real polynomials has focussed on algo- 
rithms that are quasi-linear in the degree. See, e.g., [EMOl].) 



Echoing the parallels between FEASr and ADISC> pro- 
vided by Theorems IQ1 and [T71 both FEAS R (^"1,4) GP and 

7 

ADISC> (-7-1,4) 6 P are open problems. 

As for earlier complexity lower bounds for FEASr in 
terms of sparsity, we are unaware of any. Indeed, it is not 
even known whether FEASr(Z[2;i, . . . , x n ]) is NP-hard for 
some fixed n. Also, complexity lower bounds for the vanish- 
ing of discriminants of n-variate (n + fc(n))-nomials (with k 
a slowly growing function of n) appear to be new. However, 
recent work shows that the geometry of discriminants cham- 
bers can be quite intricate already for / £^3,3+3 DRRS07 . 
Also, it was known even earlier that deciding the vanish- 
ing of sparse discriminants of univariate m-nomials (with 
m unbounded) is already NP-hard with respect to random- 
ized reductions [KS99j . Considering Theorems 11.31 and 11.71 
one may thus be inclined to conjecture that FEASi(Z[a;i]) 
is NP-hard. Curiously, over a different family of complete 
fields (thep-adic rationals), one can already prove that de- 
tecting roots for univariate m-nomials (with m unbounded) 
is NP-hard with respect to randomized reductions [IRR07] . 

2. BACKGROUND AND ANCILLARY 
RESULTS 

After recalling a basic complexity construction, we will 
present some tools for dealing with n-variate (n+1) -normals, 
and then move on to n-variate (n + fc)-nomials with k > 2. 
All proofs for the results of this section are in the 
Appendix. 

2.1 A Key Reduction 

To measure the complexity of our algorithms, let us fix 
the following definitions for input size. 

Definition 2.1. For any aGZ, we define its size, size(a), 
to be 1 + log(l + \a\). More generally, we define the size of 
a matrix U — [tH,j] G Z mxn to be J]\ . size(uij). Also, for 
any f(x) — YliLi c i x<li G Z[xi, . . . ,x n ], we define size(/) to 
be J27LA size ( c i) + size(oj)]. Finally, for F = . . . , f k ) G 
(Z[xi, . . . ,x„]) k , we define size(F) = J^ i=1 size(/i). o 

A key construction we will use later in our NP-hardness 
proofs is a refinement of an old trick for embedding Boolean 
satisfiability into real/complex satisfiability. We refer to the 
well-known 3CNFSAT problem, reviewed in the Appendix. 

Proposition 2.2. Given any 3CNFS AT instance B(X) 
with n variables and N clauses, let Wb denote 
(({1} x P£) U (Pc x {l})) 4N ~ n . Then there is an 
(8N — n) x (87V — n) polynomial system Fb with the 
following properties: 

1. B(X) is satis fiable iff Fb has a root in 
{1,2}" xWb. 

2. Fb has no more than 33JV — 4n monomial terms, 
size(Fs) = O(N), and every root of Fb in 
(Pc) 8JV_n lies in {l,2} n x Wb and is degenerate. 

Also, if we define tM^Zi, ■ ■ ■ , Zm) to be 
1 + zf +1 + ■ ■ ■ + zf { +l - (M + 1)21 • ■ • z M , then 

3. tm is nonnegative on R+ 7 , with a unique positive root 
at (1, ... , 1) that happens to be the only degenerate root 
oft M m C M . 



llfe>0,fe K,n+k, and M := [fc 1 ^] , then f(x) + 
tM(z)sJ-\* + s for r\ = n J r M and some positive S <e. 
In particular, size(/(x) + (m(z)) = 0^size(/) 1 ^ e ^ . ■ 

The seemingly mysterious polynomial tM defined above will 
be useful later when we will need to decrease the difference 
between the number of terms and variables in certain poly- 
nomials. 

2.2 Efficient Linear Algebra on Exponents 

A simple and useful change of variables is to use 
monomials in new variables. 

Definition 2.3. For any ring R, let R mx 71 denote the set 
of m x n matrices with entries in R. For any M = [mij] G 
M nxn and y = (yi, . . . , y n ), we define the formal expression 

U M ~ ivT 1 ' 1 ■ ■ ■ Vn"" 1 , uT 1 '" ■ ■ ■ yn"' n )- We call the sub- 
stitution x:=y M a monomial change of variables. Also, 
for any z := (21, . . . , z n ), we let xz := (xiZi, . . . , x n z n ). Fi- 
nally, let GL„(Z) denote the group of all matrices in Z nxn 
with determinant ±1 (the set o/unimodular matrices), o 

Proposition 2.4. (See, e.g., \LRW03l Prop. 2].) For 
any U, V G R nxn , we have the formal identity (xy) uv — 
(x u ) v (y u ) v . Also, if detU 0, then the function 
eu(x) := x u is an analytic automorphism o/R", and pre- 
serves smooth points and singular points of positive zero sets 
of analytic functions. Moreover, ifdetU>0, then eu in fact 
induces a diffeotopy on any positive zero set of an analytic 
function. Finally, U G GL„(R) implies that % X (R+) = R+ 
and that eu maps distinct open orthants of R™ to distinct 
open orthants o/R™. ■ 

Proposition 12.41 with minor variations, has been 
observed in many earlier works (see, e.g., [LRW03 ). Per- 
haps the only new ingredient is the observation on diffeotopy, 
which follows easily from the fact that GL^ (R) (the set of all 
n x n real matrices with positive determinant) is a connected 
Lie group. 

Recall that the affine span of a point set _4cR™, AfL4, 
is the set of real linear combinations ~^2 aeA c a a satisfying 

Lemma 2.5. Given any f G J- n ,m with d — 
dim Aff (Supp(/)) < n < m, we can find (using a number 
of bit operations polynomial in size(Supp(/)) ) aU GGL„(Z) 
such that g(y) :=f(y u ) £J~^ m and g vanishes in R+ (resp. 
(R*) d ,) iff f vanishes in R+ (resp. (R*)™^. In particular, 
there is an absolute constant c such that size(£7) = 
0(size(Supp(/)) c ). 

To study Z+(f) when / G J-n, n +i it wm help to have a 
much simpler canonical form. In what follows, we use # for 
set cardinality and d for the i— standard basis vector of R n . 

Lemma 2.6. For any f G J~n, n +i we can compute I G 
{0, . . . , n} within NC and 7 G R+ such that f(x) :— 7 + 
xi + ■ ■ ■ +xi — xe+i — • • • — x n satisfies: (1) either f or —f 
has exactly £+1 positive coefficients, and (2) Z+(f) and 
Z+(f) are diffeotopic. 

Corollary 2.7. Suppose f G J^ >n +i °-nd Supp(/) = 
{ai,...,a„ +1 }cR". Then 



1. f has a root in R+ <=>■ not all the coefficients of f 
have the same sign. In particular, Z+(f) is diffeotopic 
to either R+ _1 or 0. 

2. If all the coefficients of f have the same sign, then f 
has a root in (R*) n there are indices i € [n] and 
j,j £ [n + 1] with a,ij — a,i ji odd. 

2.3 Combinatorics and Topology of Certain A- 
Discriminants 

The connection between topology of discriminant comple- 
ments and computational complexity dates back to the late 
1970s, having been observed relative to (a) the membership 
problem for semi-algebraic sets |DL79| and (b) the approx- 
imation of roots of univariate polynomials Sma87 . Our 
goal here is a precise connection between FEASr and A- 
discriminant complements. (See also DRRS07 for further 
results in this direction.) 

Definition 2.8. \GKZ94\ Ch. 1 & 9-11] Given any A = 

{ai, . . . , a m } C Z n of cardinality m and ci, , c m £ C*, we 

define X7a C ft; 1-1 the ,4-discriminant 

variety to be the closure of the set of all 

[ci : ■ ■ ■ : Cm] GP™^ 1 such that f(x)=Y%L 1 CiX ai has a de- 
generate root in C n . We then define Aa €Z[ci, . . . , c m ]\{0} 
— the ^-discriminant — to be the unique (up to sign) ir- 
reducible defining polynomial of V a- Also, when Va has 
complex codimension at least 2, we set Aa to the constant 
1. For convenience, we will sometimes write A.a(/) in place 

of Aa{_Ci, . . . , Cm). O 

To prove our results, it will actually suffice to deal with a 
small subclass of ,4-discriminants. 

Definition 2.9. We call A C R n a (non-degenerate) 

circuit iff A is affinely dependent, but every proper sub- 
set of A is affinely independent. Also, we say that A is a 
degenerate circuit iff A contains a point a and a proper- 
subset B such that a£B, A \ a is affinely independent, and 
B is a non-degenerate circuit, o 

For instance, both \_i and 1 * are circuits, but 1 . . is a 
degenerate circuit. In general, for any degenerate circuit A, 
the subset B named above is always unique. 

The relevance of „4-discriminants to m-nomial zero sets 
can be summarized as follows. 

Definition 2.10. Following the notation of Definition ] 2. 8i 
we call any connected component o/P^ 1-1 \ (Va U{[xi : •■■ : 
Xm\ X\ ■ • ' Xm = 0}) a (real) „4-discriminant chamber. 
Also, given any subsets X,ycE™, we say that they are iso- 
topic (resp. diffeotopic ) iff there is a continuous (resp. dif- 
ferentiable) function H : [0, 1] x X — > R+ such that H(t, •) 
is a homeomorphism (resp. diffeomorphism) for all t€ [0, 1], 
H(0,-) is the identity on X, and H(1,X) = Y . Finally, for 
any A C R™ of cardinality m, let Ta denote the set of all 
n-variate m-nomials with support A. o 

Remark 2.11. Note that when A has cardinality m, we 
may naturally identify elements of P™~ (resp. PJJ 1- ) with 
equivalence classes determined by nonzero complex (resp. 
real) multiples of elements of J- a- o 

2 This terminology comes from matroid theory and has noth- 
ing to do with circuits from complexity theory. 



The topology of toric real zero sets is known to be constant 
on discriminant chambers (see, e.g., [GKZ941 Ch. 11, Sec. 
5A, Prop. 5.2, pg. 382]). However, we will need a refinement 
of this fact to positive zero sets, so we derive this directly for 
A in sufficiently general position — non-degenerate circuits 
in particular. 

Lemma 2.12. Following the notation above, suppose Ac 
R n is such that the minimum of any linear form on A is 
minimized at no more than n+1 points. Also letC be any A- 
discriminant chamber. Then f,g£C => Z+ (/) and Z+ (g) 
are diffeotopic. 

There is then a very compact description for when A 
is a circuit. 

Lemma 2.13. Suppose A= {ai, a„+2} CZ™ is a non- 
degenerate circuit, f is a polynomial with support A, A is 
the (n+1) x (n+2) matrix whose j— column is {1} x a 3 ■, Aj is 
the submatrix of A obtained by deleting the j— column, and 
bj :— ( — l) 3 det Aj/f3 where [3 = gcd(det Ai , . . ■ , det An+2) ■ 
Then: 

1. A(ci, . . . , Cn+2) is, up to a multiple by a nonzero mono- 

n + 2 , sbi 

mial term, Yl ( ) — 1- Also, (61, . . . , 6^+2) can be 
computed in P. 

n + 2 

2. n (sign(6 I c0c I /60 slsn(, ' lC,) ' >i = 1 for some [ci : • • • : 

c, l+ 2]ePS +1 withsign(c 1 b 1 )= ■■■ =sign(c n+2 bn+2) 
Z+ (X^"=i 2 c* 3 -"* ) contains a degenerate point In par- 
ticular, Z+(f) has at most one degenerate point. 

3. A has exactly two triangulations: one with simplices 
{Conv(.A \ {bi}) I sign(6i) > 0}, and the other with 
simplices {Conv(^4 \ {bi}) \ sign(fei) < 0}. Moreover, 
the preceding description also holds when A is a de- 
generate circuit. 

2.4 Complexity of Circuit Discriminants and 
Linear Forms in Logarithms 

Theorem 11.71 is a central tool behind the upper bounds 
and lower bounds of Theorem 11.31 and is precisely where 
diophantine approximation enters our scenery. To wit, the 
proof of Assertion (1) of Theorem 11.71 makes use of the fol- 
lowing powerful result. 

Nesterenko-Matveev Theorem. \NesQ3[ Thm. 2.1, Pg. 
55] For any integers c 1 ,ai, . . . ,cn,olm with ai > 2 for all 
i, define A(c, a) := cilog(ai) + ■■■ + Cjvlog(ajv). Then 

A(c, a) 7^0 log j A ^ a j j is bounded above by 

N 

2.9(iV + 2) 9/2 (2e) 2iv+6 (2 + logmax J | Cj |) H log|cy|. ■ 

j=i 

Assertion (1) of Theorem 11.71 will follow easily from the 
two algorithms we state below, once we prove their correct- 
ness and verify their efficiency. However, we will first need 
to recall the concept of a gcd-free basis. In essence, a gcd- 
free basis is nearly as powerful as factorization into primes, 
but is far easier to compute. 

Definition 2.14. fSSM , Sec. 84] For any subset {ai, . . . ,ajv} 
CN, a gcd-free basis is a pair of sets ({74}^, {ey }( 1 ,j)g[jv]x[r ) ]) 
such that (1) gcd(7^,7j) = 1 for all i ^ j, and (2) ai 
1 1' 7 e ^ for all i. o 



Algorithm 2.15. 
Input: Integers ai, iti, vi, . . . , on,I3n, un, vn- 
Output: A true declaration as to whether a" 1 • • • a^ N 



Description: 

0- //nf =1 (sign«, 



\ui mod 2 



output 'They are not equal." and STOP. 



then 



1. Replace the cti and (3i by their absolute values and then 
construct, via Theorem \4-9\ of the Appendix, a gcd- 
free basis ({7i}j=i, { e *j}(ij)e[2iv] x [r,]) for ai,...,a N , 

2- If 121 

=1 djUi — X]j=jv+i e ij v i f or M j £ M then output 
'They are equal. " and STOP. 

3. Output 'They are not equal . " 

Algorithm 2.16. 
Input: Positive integers a\,Ui,. 
. . . , Pn,vn with ai, /3i>2 for all i. 
Output: The sign of a" 1 ■ ■ ■ a^ 1 - 
Description: 

0. Check via Algorithm \MB whether a" 1 ■ • ■ a\f = /3" 1 
If so, output "They are equal. " and STOP. 



a M ,u M and fh,Vi, 



• B VN 



, vn}, and 
(1 + log CO 



1. Let U :=max{ui, . . . ,um,Vi, 

M \ f N 

n log km n log i a 

2. For all i £ [M] (resp. i € [N]), let Ai (resp. Bi) be a 
rational number agreeing with logo^ (resp. \ogf3i) in 
its first 2 + E + log 2 M (resp. 2 + E + log 2 N ) leading 
bits^ 

/ M \ / N 

3. Output the sign of I UiAi J — I E v iBi 



and STOP. 



Lemma 2.17. Alaorithms \2.15\ and \2.16[ are both correct. 
Moreover, following the preceding notation, Alaorithms \2.15\ 
and \2.16\ run within a number of bit operations asymptoti- 
cally linear in, respectively, 



E (log(w») log(aO + log(t»i) log(A)) 2 



and 



(M + N)(30) M+N L(lagU) fl LQaeM) 



n i(iog(ft)) 

(=1 



where L(x):=x log (x) loglog(a;). 



2.5 Positive Feasibility for Circuits 

For a real polynomial supported on a non-degenerate cir- 
cuit, there are just two ways it can fail to have a positive 
root: a simple way and a subtle way. This is summarized 
below. Recall that the Newton polytope of / is sim- 
ply Newt(/) := Conv(Supp(/)), where Conv(S') denotes the 
convex hull (smallest convex set) containing S. 

Theorem 2.18. Suppose f(x) = c^" 4 S ^,„ +2 , 

Supp(/) is a non-degenerate circuit, and b is the vector from 
Lemma \2.13\ Then Z+(f) is empty iff one of the following 
conditions holds: 



3 For definite ness, let us use Arithmetic-Geometric Mean It- 
eration as in Ber 03 to find these approximations. 



1. All the Ci have the same sign. 

2. Newt(/) is an n-simplex and, assuming ay is the unique 
element of A lying in the interior of Newt(/), 
we have — sign(cj/) =sign(ci) for all i=£j' and 



ns 2 (ssgnfc 



>i. 



Positive feasibility for polynomials supported on degen- 
erate circuits can then essentially be reduced to the non- 
degenerate circuit case in some lower dimension. An addi- 
tional twist arises from the fact that the zero sets of polyno- 
mials supported on degenerate circuits are, up to a monomial 
change of variables, the graphs of polynomials supported on 
non-degenerate circuits. 

Theorem 2.19. Suppose f(x) = J27=i CiX,ai £?n,n+2 has 
support AcV. n that is a degenerate circuit with non- degenerate 
subcircuit B = {ai, . . . , ay}, and b is the vector defined in 
Lemma \2.13\ (ignoring the non- degeneracy assumption for 
A). Then, when not all the coefficients of f have the same 
sign, Z+(f) is empty iff both the following conditions hold: 
a. Conv(S) is a (j' — 2)-simplex and, permuting indices so 

that ay is the unique element of B lying in the relative 

interior of Conv(S), we have — sign(cy) = sign(ci) for 

all ij^j'. 

6- UU (sfen&<*)£ ) >1. 

3. THE PROOFS OF OUR MAIN RESULTS: 
THEOREMS O AND O 

We go in increasing order of proof length. 

3.1 Proving Theorem 11.71 

Assertion (1): First note that any input / must have sup- 
port A — {ai, . . . ,a n +2} equal to either a degenerate cir- 
cuit or a non-degenerate circuit. Recalling Assertion (1) of 
Lemma [2.13l observe then that the vector b:= (bi, . . . , bn+2) 
has a zero coordinate iff A is a degenerate circuit, and 6 
can be computed in time polynomial in size(,4). If A is a 
degenerate circuit then (following easily from the definition) 
A^i must be identically 1, thus leaving Assertion (1) of our 
present theorem trivially true. So let us assume henceforth 
that A is a non-degenerate circuit, and that Cj is the coeffi- 
cient of x a i in / for all j. 

Via Assertion (1) of Lemma 12.131 once again, Assertion 
(1) of Theorem 11.71 follows routinely from the complexity 
bounds from Lemma 12.171 In particular, the latter lemma 
tells us that the bit complexity of ADISC = , for input coef- 
ficients (ci, . . . , Cn+2), is polynomial in J27=i logte&O (fol- 
lowing the notation of Lemma I2.13P ; and the same is true 
for ADISC> provided n is fixed. The classical Hadamard 
inequality [Mig92 then tells us that size(fci) is 
0(n log(n maxj^jajfc})). So the complexity of ADISC = 
is indeed polynomial in size(/); and the same holds for 
ADISC> when n is fixed. ■ 

Assertion (2): We will construct an explicit reduction of 
3CNFSAT to ADISC= ( (J Tl in+n e ) . In particular, to 

VnGM / 

any 3CNFSAT instance B{X) with N clauses and n vari- 
ables, let us first consider Fb = (/i, . . . , fsN-n) — the asso- 
ciated (8JV — n) x (8N — n) polynomial system as detailed 
in Definition 14.11 of the Appendix and Proposition 12.21 of 
Section O 



Let us then set M := | (max{0, 17 N - 2n} + 2) 1/£ | and 
define the single polynomial fs to be 

fl + Al/2 + ■ ■ ■ + XsN-n-lf&N-n + AgJV -nt M (Zl, . . . , Zm) ■ 

Letting A be the support of fn, it is then easily checked 
(from Definition 14.11 of the Appendix and Proposition I2,2|l 
that A is affinely independent and /a is in J-i 6N _ 2n+M N > 
for some N' < 33N - An + M + 2. 

By the Cayley Trick [GKZ941 Prop. 1.7, pp. 274] we 
then obtain that Ayi(/s)=0 iff 

(*) Fb has a degenerate root in (Pc) 2JV_n and tu has 
a degenerate root in (C*) AI . 
(Since Newt(tjvf) is a simplex, it is easily checked that tu 
has no complex degenerate roots at infinity.) By Proposition 
12.21 the degenerate roots of Fb are exactly {1, 2} n x Wb, 
and Im has a unique degenerate root by construction. So 
(*) holds iff B(X) has a satisfying assignment. We have 
thus reduced 3CNFSAT to detecting the vanishing of a 
particular „4-discriminant. 

To conclude, observe that the number of terms of /s is 
only slightly larger than its number of variables, thanks 
to Proposition 12.21 In particular, size(/s) = 0(size(B) 1 ^ £ ) 
and f B G UngN-^n.n+n* for some S e (°i e ]- Clearly then, 
ADISC = (U neN J r i,„+„0 G P => P = NP, thus proving 
our first desired NP-hardness lower bound. 

The NP-hardness of ADISC> ([} CM T* , s ) then fol- 

lows immediately since ADISC> ^U„gM -Fn,n+n s ) is a re ~ 
finement of ADISC= ( M _ W .F* , s ). U 

3.2 Proving Theorem [L3] 

Assertion (2): We will give an explicit reduction of 3CNFSAT 
to FEAS+ (UngN 3~n,n+n E ) • Attaining such a reduction will 
require little effort, thanks to our earlier reduction used to 
prove Assertion (2) of Theorem 11.71 

In particular, for any 3CNFSAT instance B with N 
clauses and n variables, let us recall the system 
Fb = . . . , fsN-n) from Definition 14.11 of the Appendix 
and Proposition 12.21 Let us then define M to be 

(max{0, 42AT - n} + 2) 1/£ "| and define g B (x,z) to be 

fi(x) + ■ ■ ■ + fi N (x) + t M {zi, ■ ■ ■ , z M )- It is then easily 
checked that f B G ^„+m.n' for some N ' < + A/ + 2. 
Moreover, B has a satisfying assignment iff qb has a posi- 
tive root. (Indeed, any root of qb clearly lies in {1,2}™ x 
{1} M .) We have thus reduced 3CNFSAT to a special case 
of FEAS+. 

Now observe that the number of terms of is only slightly 
larger than its number of variables, thanks to Proposition 

12.21 In particular, size(g_B) = 0^size(£>) 1//e j and gs is in 
U„ s n ^,„ +n i Ior some 3 G (Ml- Clearly then, 
FEAS+(U neN ^,„+„ e ) 6 P => P = NP, thus proving 
one of our desired NP-hardness lower bounds. 

To conclude, we now need to prove the NP-hardness of 
FEASi(U„ €M f,t,„ + „e). This we do by employing our pre- 
ceding argument almost verbatim. The only difference is 
that we instead use the polynomial hs{x, z) :=/i (x) + ■ ■ ■ + 
IIn{x) + t M (zi, Zm), and observe that t M (zf, . . . , z 2 M ) 
is nonnegative on all of R n . So we are done. ■ 

Assertion (0): Our topological assertion follows immedi- 
ately from Lemma [2. 61 and Corollary 12.71 

To obtain our algorithmic assertions, simply note that by 



Assertion (1) of Corollary 12.71 detecting positive roots for 
/ reduces to checking whether all the coefficients have the 
same sign. This can clearly be done by n sign evaluations 
and n — 1 comparisons, doable in logarithmic parallel time. 
So the inclusion involving FEAS+ is proved. 

Let us now show that we can detect roots in (R*) n within 
NC 1 : Employing our algorithm from the last paragraph, 
we can clearly assume the signs of the coefficients of / are 
all identical (for otherwise, we would have detected a root in 
M™ and finished). So then, by Assertion (2) of Corollary |2.7l 
we can simply do a parity check (trivially doable in NC 1 ) 
of the entries of [012 — ai, ... , a n +i — ai]. 

To conclude, we simply observe that our algorithm for 
detecting roots in (R*) n trivially extends to root detection 
in R n : Any root of / in R n must lie in some coordinate 
subspace L of minimal positive dimension. So, on L, the 
honest n-variate (n + l)-nomial / will restrict to an /' 6 
JF*, n / +1 with n' < n and support a subset of the columns 
of a submatrix of A. So then, we must check whether (a) 
all the coefficients of /' have the same sign or (if not), (b) 
a submatrix of [a 2 — ai, . . . ,a n +i — a\] has an odd entry. 
In other words, / has a root in R n <^=> / has a root in 
(R*) n U {O}. Since checking whether / vanishes at O is 
the same as checking whether / is missing a constant term, 
checking for roots in R n is thus also in NC 1 . ■ 

Remark 3.1. Note that checking whether a given f 6 
■F„, n +i lies in J~n, n +i can be done within NC 2 / one sim- 
ply finds d = dim Supp(/) in NC 2 by computing the rank 
of the matrix whose columns are a 2 — ai, . . . ,a m — ai (via 
the parallel algorithm of Csanky \Csa70j ), and then checks 
whether d — n. o 

Assertion (1): The algorithm we use to prove FEAS + (J r * ri _|_ 2 ) 
G P for fixed n is described just below. Note also that once 
we have FEAS+(jF^ in+2 ) G P for fixed n, it easily follows 
that FEASr(JT* n+2 ) 6 P: The polynomial obtained from 
an / G ^Fn,n+2 by setting any non-empty subset of its vari- 
ables to clearly lies in JF*, n , +2 for some n' <n (modulo a 
permutation of variables). Thus, since we can apply changes 
of variables like Xi 1— » — Xi in P, and since there are exactly 
3 n sequences of the form (ei, . . . , e n ) with Si G {0, ±1} for all 
i, it thus clearly suffices to show that FEAS+(J r ^ n+2 ) GP 
for fixed n. 

We thus need only prove correctness, and a suitable com- 
plexity bound, for the following algorithm: 

Algorithm 3.2. 
Input: A coefficient vector c := (ci, . . . , c n + 2 ) and a 
(possibly degenerate) circuit A = {ai, . . . , a n + 2 } of cardinal- 
ity n + 2. 

Output: A true declaration as to whether Z+(f) is empty 

or not, where f(x):—'^2"^CiX ai . 

Description: 

1. If all the d have the same sign then output "Z+(/) = 0" 
and STOP. 

2. Let b = (61, ... , 6 n + 2 ) G be the vector obtained by 
applying Lemma \2.1'A to A. If b or —b has a unique 
negative coordinate bji , and Cj> is the unique negative 
coordinate of c or —c, then do the following: 

(a) Replace b by — sign(6j/)6, replace c by — sign(c 3 -/)c, 
and then reorder b, c, and A by the same 
permutation so that bj/ < and [bi>0 iff i< j']. 



(b) Ifj'<n + 2 and 

v' ii: <■;■ nri 1 ^ 

then output "Z + (f) = <&" and STOP. 

(c) Decide via Alaorithm \2.16\ whether 

bp' n;',v '-/'nr 

If so, output "Z+(f) = Q" and STOP. 

Output "Z+(f) is non-empty! " and STOP. 

The correctness of Algorithm 13.21 follows directly from The- 
orems [2TT9] an d [2TT8J In particular, note that bi is simply the 
signed volume of Conv(.A\{ai}). So the geometric interpre- 
tation b or —b having a unique negative coordinate is that 
the convex hull of the unique non-degenerate subcircuit of A 
is a simplex, with ay lying in its relative interior. Similarly, 
the geometric interpretation of j' < n+2 is that A is a degen- 
erate circuit. Finally, the product comparisons from Steps 
(b) and (c) simply decide the product inequalities stated in 
Theorem |2~T81 and Theorem ETT91 

So now we need only bound complexity, and this follows 
immediately from Lemma [2.17l (assuming we use Algorithm 
EH] for Step (b)). ■ 

It is worth noting that we need to compute the sign of a 
linear combination of logarithms only when the unique non- 
degenerate subcircuit B of A is a simplex, and all "vertex" 
coefficients have sign opposite from the "internal" coefficient. 
Also, just as in Remark 13.11 checking whether a given / G 
J~n,n+2 lies in T^n+s can be done within NC 2 by computing 
d = dimSupp(/) efficiently. Moreover, from our preceding 
proof, we see that deciding whether a circuit is degenerate 
(and extracting B from A when A is degenerate) can be 
done in NC 2 as well, since we can set (3= 1 if we only want 
the signs of (61, . . . , b n+2 ). 
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Appendix: Complexity, Viro Diagrams, and Post- 
poned Proofs 

4.1 Complexity Classes and 3CNfsat 

A complete and rigourous description of the complexity 
classes we used can be found in [Pap95| . So for the conve- 
nience of the reader, we briefly review the following defini- 
tions: 

NC The family of functions computable by Boolean 
circuitifl with size polynomial in the input size and 
depth 0(log l InputSize). 

P The family of decision problems that can be done within 
time polynomial in the input size. 

NP The family of decision problems where a "Yes" answer 
can be verified within time polynomial in the input 
size. 

PSPACE The family of decision problems solvable within 
polynomial-time, provided a number of processors 
exponential in the input size is allowed. 

Definition 4.1. Recall that 3CNFSAT is the problem 
of deciding whether an n-variate Boolean formula of the form 
B(X) — Ci(X) A ■ ■ ■ A Cn{X) has a satisfying assignment, 
where each clause Ci is of one of the following forms: 
X l \fX j yX k , XjV Xj\/^X k , 
Xi V -iXj V -^X k , -^Xi V -lXj v -*x h , 

i,j,k £ [n] are pair-wise distinct, \^~\ < N < 8 > an< ^ 

a satisfying assigment consists of an assignment of values 
from {True, False} to the variables Xi, . . . ,X n yielding the 
equality B(X) = True 'GJ79I. We then define size(B) := 

4 This is the one time we will mention circuits in the sense 
of complexity theory: Everywhere else in this paper, our 
circuits will be combinatorial objects as in Definition 



3N, a(xi,x 2 ,x 3 ) := (xi — 2)(x 2 — 2)(x 3 - 2), and b(x 1 ) := 
(xi — l)(xi — 2). Finally, to any 3CNFSAT clause Ci as 
above, we associate a 4 x 3 polynomial system Hc e as fol- 
lows: we respectively map clauses of the form Xi V Xj V Xk, 
X z V Xj V X z V -iX,- V -iX fc) -.X; V -iXj V ^X fe to 

quadruples of the form 

(a(xi, Xj,Xk), b(xi), b(xj), b(xk)), 
(a(xi,Xj,3- Xk),b(xt),b(xj),b(xk)), 
(a(xi,3- Xj,3- Xk),b(xi).b(xj),b(x k )), 
(o(3 - a;,, 3- Sj,3 - Xk),b(Xi),b(xj),b(xk)); 
and we associate to the 3CNFSAT instance B(X) a 
(8N — n) x (8N — n) polynomial system with integral co- 
efficients, Fb, defined to be 

(%,. • ■ , H Cn , (Wl-l)(«l — 1), . . . , (U4JV- n — l)(v4JV-n— 1)). 

In particular, assigning True (resp. False,) to Xi will corre- 
spond to setting Xi = 2 (resp. Xi = l). o 

Note that Fb has a natural and well-defined zero set in 
(Pc) 8JV_n since its Newton polytopes are either axes-parallel 
line segments or 3-cubes, and we can multihomogenize with 
8N — n extra variables. 

Proof of Proposition l2.2t Assertions (1) and (2) of Propo- 
sition |2~21 are elementary. In particular, the last 4iV — n poly- 
nomials of Fb simply ensure that Fb has enough variables 
so that it is square. Assertion (3) follows easily from the 
classical Arithmetic-Geometric Inequality IHLP881 Sec. 2.5, 
pp. 16-18]. Assertion (4) follows easily upon observing the 

inequalities k < = M E < (n + Alf and the fact that 

Newt (4m) is M- dimensional. ■ 

4.2 Digression on Viro Diagrams 

Let us recall an elegant result of Oleg Viro on the clas- 
sification of certain real algebraic hypersurfaces. In what 
follows, we liberally paraphrase from Proposition 5.2 and 
Theorem 5.6 of [GKZ94I Ch. 5, pp. 378-393]. 

Definition 4.2. Given any finite point set A C R™, let 
us call any function ui : A — ► R a lifting, denote by 
7r : R n+1 — > W 1 the natural projection which forgets the 
last coordinate, and let A := {(a,u(a)) | a £ A}. We then 
say that the polyhedral subdivision E w of A defined by 
{tt(Q) I Q a lower facet of Conv A of dimension dimAfL4} 
is induced by the lifting u), and we call E u a triangula- 
tion induced by a lifting iff every cell of E u is a simplex. 
Finally, given any f(x) = ^2 a ^^ c aX a £ Z[x±, . . . , x„], we 
define / w , e (x) := Y2 a eA CaS^^x 0, to be the toric perturba- 
tion of / (corresponding to the lifting u). o 

Definition 4.3. Following the notation above, suppose 
dimAff.4 = n and A is equipped with a triangulation E in- 
duced by a lifting and a function s : A — > {±} which we 
will call a distribution of signs for A. We then locally 
define a piece-wise linear manifold — the Viro diagram 
Va(E,s) — in the following local manner: For any n-cell 
C £ E, let Lc be the convex hull of the set of midpoints of 
edges of C with vertices of opposite sign, and then define 
V.a(E,s) := (J Lc- When A — Supp(/) and s is the 

C an n-cell 

corresponding sequence of coefficient signs, then we also call 
V(f) := V.a(E, s) the Viro diagram of /. o 



5 A lower facet is simply a facet which has an inner normal 
with positive last coordinate. 



Example 4.4. The following figure illustrates 6 circuits 
of cardinality 4, each equipped with a triangulation induced 
by a lifting, and a distribution of signs. The corresponding 
(possibly empty) Viro diagrams are drawn in thicker lines, o 




Vino's Theorem. Suppose f(x) = J2 a eA CaXa is in 
R[xi, . . . , x n ] with Supp(/) = A and dimAff.4 = n, ui is 
any lifting of A, and define s/(a) = sign(c a ) for all a £ A. 
Then for any sufficiently small e>0, Z+(f UJ}£ ) is isotopic to 
Va(^ui , s f)\dConvA. In particular, VU(E w ,s/) is a disjoint 
finite union of piece-wise linear manifolds, each possibly hav- 
ing a non-empty boundary. M 

Lemma 4.5. Suppose A is a circuit, E is a triangulation 
of A, n = dimAff.4, and s is any distribution of signs on A. 
Then Z+(f) smooth => Z+(f) is isotopic to V.a(E,s/) \ 
9Conv^4 for some E. 

Proof of Lemma l4.5l By Lemma[27l3]it easily follows that 
A has at most 2 discriminant chambers in R" +2 , and each 
such chamber contains a unique toric perturbation. Since 
the topology of Z+(f) is constant on any discriminant cham- 
ber containing / by Lemma 12.121 we are done. ■ 

Another important consequence of Viro's Theorem deals 
with how "roots at infinity" sometimes imply the existence 
of non-compact connected components for a positive zero 
set. 

Lemma 4.6. Assume f has support such that every facet 
o/Newt(/) is a simplex and coefficients such that V{f) in- 
tersects 9Newt(/) for every underlying triangulation. Then 
Z+(f) has a non-compact connected component. ■ 

Employing our current notation, Lemma [4 . 6 1 follows directly 
from Lemma 15 of |LRW03| . 

4.3 More Postponed Proofs 

Proof of Lemma 12.51 Lemma 12.51 follows immediately 
from the following well-known factorization for integer ma- 
trices and its recent complexity bounds. 

Definition 4.7. JIM3L[StbWf Given anu M £Z mx " . the 
Hermite factorization of M is an identity of the form 
UM = H where U £ GL,„ and H = [hi 3 ] £ Z nxn is non- 
negative and upper triangular, with all off-diagonal entries 
smaller than the positive diagonal entry in the same col- 
umn. Finally, the Smith factorization of M is an iden- 
tity of the form UMV — S where U £ GL m , V £ GL n , and 
S= [sij] £ Z mx " is diagonal, with Si t i\si+i,\+i for alii, o 

Lemma 4.8. UhM \Sto9$ For any M = [m, ,] £ Z" x " . the 
Hermite and Smith factorizations of M exist uniquely, and 
can be computed within 0(n 4 log 3 (nmaxi,j bit opera- 

tions. Furthermore, in the notation of Definition ^. 7[ the en- 
tries of U, V, S, and H all have bit size 
0(n 3 log 2 (2n + maxij \mij\)). ■ 



Proof of Lemma 12.61 Computing our desired canonical 
form / boils down to reordering monomials, performing a 
monomial change of variables, and a rescaling. 

First, let us replace / by sign(ci)/ and reorder the terms 
of / supported on {<J2, • • • , a„ + i} so that 02, . . . , <v (resp. 
C(/ + i, . . . , Cn+i) are positive (resp. negative), for some unique 
£' £ {1, • • • , n + 1}. We then form the n x n matrix B whose 
whose i— column is a i+1 — oi, for i £ {1, . . . , n}. If neces- 
sary, let us also swap the terms supported on {02,03} (or 
{a n ,a n +i}), before defining B, so that det B > and the 
sign condition defining £' is still preserved. (If n — 1 then 
we can simply reorder terms so that ai < a-2 and define B 
accordingly.) 

Letting fix) := 1 — — we then obtain by Proposition 

x a l 

12.41 that Z+(f) and Z+(jf^j are diffeotopic. Moreover, we 
clearly have that 

f(x) :—Ci + C2X\ + - ■ ■+c n +iXn (remember we have permuted 
the terms of /, and thus the Ci as well) where ci >0. So we 

can now define f(x) to be f(j^, ••■ j [^^J J 1 define 7 to be 
the constant term of /, and set £:=£' — 1. It is then easily 
verified that the coefficients of xi, . . . , xi (resp. . . . , x n ) 

are all 1 (resp. —1). 

Since / was defined by rescaling the variables of /, and 
such a scaling of variables can easily be parametrized so as to 
induce a diffeotopy from Z+ (/) to Z+ \jf\ , Assertion (2) fol- 
lows. Assertion (1) holds trivially since we never altered the 
difference between the number of positive and negative coef- 
ficients. That £ is computable in NC 1 follows from Lemma 
14.81 we can compute £ simply by sorting, which clearly lies 
in NC 1 . ■ 

Proof of Corollary [2?fl 

Assertion (1): Employing the canonical form f(x) = 7 + 
xi + ■ ■ ■ + xe — xe+i — ■ ■ ■ — x n of Lemma 12.61 (with 7 > 
and £ £ {0, . . . , n}, by construction), the desired equivalence 
will follow upon proving that / has a root in K™ •<=>■ £<n. 
The latter equivalence is trivially true. By Lemma |2 . 6 1 once 
more, the statement on diffeotopy type can be reduced to 
the special case of /, which is also immediate. ■ 
Assertion (2): Dividing by a suitable monomial term, we 
can clearly assume that a\ = O and all the coefficients of 
/ are positive. So it suffices to prove that / has a root in 
(K*) n there are indices i £ [n] and j £ {2, . . . , n + 1} 

with dij odd. Writing f(x) = a + c 2 x a2 H + c n+ ix a ' n + 1 , 

let us now prove the last equivalence. 

(=>): To prove the contrapositive, assume every a%,j is 
even. Then the sign of / is positive on all of (R*) n and / 
thus has no roots in (K*) n . 

(•<=): Reordering terms and variables, we can clearly 
assume that 01,1 is odd. Letting g(x±) := f(xi, 1, . . . , 1), 
note that g must be of the form g(x\) =Cj + c^x" 1 ^ , where 

ci,C2 >0. Since g has the well-defined real root — 



then clearly has the root 
(R*) n . ■ 



1 I which lies in 



Proof of Lemma l2.12t First, recall that subanalytic sets 
are those sets defined by projections of feasible sets of sys- 
tems of analytic inequalities. In particular, C is a subanalytic 
set, and C is thus path connected since C admits a decom- 



position into connected cells. The existence of such a cell 
decomposition follows immediately from the o-minimality 
of subanalytic sets vdD98, vdD86 . Moreover, C is path 
connected via differentiable paths by the classical density 
of C u functions among C 1 functions (see, e.g., [Hir941 Ch. 
2])- 

So let m = #A and let <j> := [fa : ■ ■ ■ : </> m ] : [0, 1] — > C 
be any differentiable path connecting / and g. Also let T 
be the positive part of the real toric variety corresponding 



to Conv(A), 1a :={(c,a;)€l 



x T I E^i^ ai =0} the 



underlying (real) incidence manifold, and let ty denote the 
natural projection mapping P^ 1 " 1 x T — ► P^ 1-1 . Note then 
that cj> induces an embedded smooth compact submanifold 
M C I a, consisting of all those (c, x) with c = fat) and 
CiX ai = for some t. In particular, we see that M is 
fibered over [0, 1] and that Tp = fa 1 o ty is a Morse function 
on M with no critical points in [0, 1]. More to the point, we 
obtain a natural flow on M inducing a diffeotopy 8 between 
the zero sets of / and g in T |Hir94l Thm. 2.2, pg. 153]. 

To conclude, we simply observe that the intersection of 
, i/' _1 (t) with toric infinity is smooth for all t £ [0, 1] (by our 
assumption on the facets of Conv(^l)) and thus 8 restricts 
to a diffeotopy between Z+(f) and Z+(g). ■ 



Proof of Lemma l2.13t With the exception of the assertion 
on complexity, Lemma 12.131 follows directly from [GKZ94 
Prop. 1.8, Pg. 274], [GKZ94I Prop. 1.2, pg. 217], and the 
discussion following up to the end of Section B on page 218 of 
GKZ94]. In particular, the factor f3 takes into account that 
A may not affinely generate Z n , but is always the integral 
affine image of an A' that is. So the sign condition arises 
simply from a binomial system (with odd determinant) that 
C must satisfy. 

The assertion on the complexity of computing {pi, ... , b n +2) 
follows immediately from Csanky's famous parallel algorithm 
for the determinant |Csa76] . combined with Lemma 14.81 
Indeed, were it not for the gcd computation for j3, we could 
instead assert an NC 2 complexity bound. ■ 

Proof of Lemma I2.17t We first recall the following theo- 



Theorem 4.9. \BS90, Thm . 4.8.7, Sec. 4.8] Following 
the notation of Definition \2.14\ there is a gcd-free basis for 
{ai, . . . , ajv}, withrj, size(7;), and size(eij) each polynomial 
in y\_,| size(a^), for all i and j. Moreover, one can always 

find such a gcd-free basis using just Q^ElLi size(a?i)^ 
bit operations. ■ 



Returning to the proof of Lemma 12.171 note then that 
Algorithm 12.151 is correct and runs in the time stated by 
Theorem 14.91 and the naive complexity bounds for integer 
multiplication. 

To prove the remaining half of our lemma, observe that 



the sign of a" 1 ■ ■ ■ a\ 



Pj 1 ■ ■ ■ is the same as the 



sign of S := (X^WilogaiJ - (j2^ =1 Vi log fit J . Clearly 

then, |S- [(Etei«t>4*) - (££i«i#i)]| < E/2 by the 
Nesterenko-Matveev Theorem. So Step (3) of Algorithm 
12. 161 indeed computes the sign of S and we thus obtain cor- 
rectness. 

To see that Algorithm 12.161 runs within the time stated, 
first note that the algorithm computes M (resp. N) approxi- 
mations of logs of positive integers, each of size 0(maxi log |ai|) 



(resp. 0(maxj log |/3i|)), correct in their first 0(E + log M) 
(resp. 0(J5 + log N)) leading bits. Employing the explicit bit 
complexity estimates for fast multiplication of [BS96 , Table 
3.1, pg. 43], it is easily checked that Bernstein's method 
quoted above uses 0(b log 2 (fe) log log b) bit operations. So, 
via our chosen method for approximating logarithms [Ber03] , 
we see that the complexity of Algorithm 12.161 is 
O (M(E + log M) log 2 (E + log M) log log(£ + log M) 

+N(E + log N) log 2 (E + log N) log log(_B + log N)) 
= 0((M + N)E log 2 (E)loglogE). Upon observing that 

M + N = o((Efeiloglog|«i|) + (Eti lo g lo glAl)), 
(2e) 2 <30, and abclog(abc) < alog(a)61og(6)clog(c) for a, b, c 
sufficiently large, our final asserted complexity bound follows 
easily. ■ 

Proof of Theorem 12.181 : First note that Condition (1) 
implies that / maintains the same (non-zero) sign through- 
out R+. So Condition (1) trivially implies that Z+(f) = 0, 
and we may assume henceforth that not all the coefficients 
of / have the same sign. 

Now, if Newt(/) is not a simplex, then every point of A 
is a vertex of Newt(/) and thus, independent of the trian- 
gulation, any Viro diagram for A must be non-empty. Since 
there are only two discriminant chambers (by Lemma [2.13p . 
Lemma [4. 51 thus implies that Z+(f) must be non-empty, as- 
suming Aa(/) 7^0. Lemma 12.131 tells us that Z+(f) must 
be non-empty if A^(/) =0. So we may assume henceforth 
that Newt(/) is a simplex and that A^(/)^0. 

Continuing our focus on Condition (2), note that if the 
sign equalities from Condition (2) fail, then there must exist 
coefficients a and of opposite sign such that ai and 
vertices of Newt(/). So, again, independent of the trian- 
gulation, any Viro diagram for A must be non-empty and 
thus (just as in the preceding paragraph) Z+ (/) must again 
be non-empty. So we may assume henceforth that the sign 
equalities from Condition (2) hold. 

At this point, it is clear that we need only show that (un- 
der our current assumptions) Z+ (/) = the discrimi- 
nant inequality from Condition (2) holds. Toward this end, 
observe that the lifting that assigns ay i— > —1 and a% i— ► for 
all i^j' induces the unique triangulation of A consisting of 
a single simplex. In particular, the underlying Viro diagram 
is empty, due to the sign equalities. So by Viro's Theorem 
and Lemma 14.51 Z+(f) is empty for |a,-/| sufficiently small. 
By Lemmata 12.121 and 12.131 this topology persists for \cji\ 
just small enough to enforce the discriminant sign stated in 
Condition (2), so we are done. ■ 

Proof of Theorem 12. 191 First, we observe that via an ar- 
gument almost identical to the proof of Lemma 12.61 we can 
find a monomial change of variables (and multiply by a suit- 
able monomial term) so that J(x) :=x v f(x M ) is of the form 

Cl + C2S" 1 H h Cji-iXjl ~ 2 2 + Cj'X a + c^+i^/l" 1 H h C n+ 2Xl n , 

where ui,...,u„eN and a€W ~ 2 x {0} n ~ j +2 . In particu- 
lar, defining fs(x) =ci + C22:" 1 + ■ ■ ■ + Cji -la^./J^ 2 + Cjix a , 
it is clear that Z+(f) is nothing more than the intersection 
of the graph of Jb (which is analytic on R^_ -2 ) with an or- 
thant. So Z+(f) is smooth and thus, by Proposition 12.41 we 
obtain that Z+(f) is diffeotopic to Z+(J). It thus suffices 
to prove our lemma for /. Note also that the conditions on 
Z+ (/) allegedly characterizing Z+ (/) = are preserved un- 
der monomial multiples and monomial changes of variables. 



Observe now that if Z+ {Jb) has > 2 points then Z+ {Jb) 
must contain a smooth point £, thanks to Lemma 12.131 So 
then, by the implicit function theorem, /e must attain both 
positive and negative values in a neighborhood of £ £ W + . 

Since the range of J B > := Cji +1 x .fs^ + ••• + c n+ 2X™ n (over 

R" - -' +2 ) must contain the positive ray or the negative ray, 
and / = Jb — (—Jb 1 )-, we thus obtain that Z+(J) is non- 
empty. So we may assume henceforth that #Z+(Jb) < 1. 

To further simplify matters, observe that if all the coef- 
ficients of Jb have the same sign then, by assumption, wc 
must have a coefficient of /g/ differing in sign from ci. So 
then, /b(R+~ 2 ) has constant sign and — /b'(R' 1_: ' +2 ) con- 
tains a ray of the same sign. Therefore Z+ (J) is non-empty. 
We may therefore also assume that not all the coefficients 
of /b have the same sign. 

We can therefore conclude by studying the two cases 
Z+(/b) = and #Z+ (/b) = 1, combined with all our 
assumptions so far. 

(£+(/e)=0): By Theorem [2T81 applied to Jb, we see 
that Conditions (a) and (b) must hold (and the inequality in 
(b) strictly so) , with the possible exception of the equalities 
— sign(cj') = sign(ci) for all i > f . So we need only show 
that Z+(/) = <^=> — sign(cj/) = sign(ci) for all i>j'- By 
an argument almost identical to the last paragraph, we can 
obtain that sign(c 3 -/) =sign(cj) for some i>f •<=>■ /b(R+) 
and — /b'QRIl - "' +2 ) intersect. So the only way we can have 
Z+(J) =0 is for the equalities — sign(Cj/) = sign(ci) to hold 
for all i>j'. 

(#Z+(/b) =1): Clearly, the sole point of Z + {Jb) must 
be singular, and thus Aa(J) ~ 0, which in turn enforces 
Condition (b) (with equality). 

Now, if Conv(B) is not a simplex, then every one of its 
faces is a simplex, since B is a non-degenerate circuit. More- 
over, all Viro diagrams of /b (regardless of triangulation) 
are non-empty and intersect dConv(B). So, by Lemma 14.61 
Z+(/b) must contain an unbounded connected component 
and thus #Z+(/b) >2. So we may assume that Conv(S) is 
a simplex. 

At this point, we need only prove that the sign equalities 
of Condition (a) must hold. Toward this end, note that 
if sign(cj/) = sign(ci) for some i £ {1, . . . , j' — 1} then all 
Viro diagrams of /b (regardless of triangulation) are non- 
empty and intersect dConv(B). So, by Lemma 14.61 again. 
Z+(Jb) must contain an unbounded connected component 
forcing #Z+(/b) > 2 again. We can therefore assume that 
-sign(c J /) = sign(c i ) for all ie{l, .. — 1}. 

To conclude, we will show that — sign(cj/)/B attains only 

nonnegative values on R^~ 2 . (This will enforce our final 

desired sign equalities sign(c J /) =sign(ci) for all i>j' 

— simply by comparing the range of /b and /b' just as 
before.) First note that /b(R+~ 2 ^ is unaffected by in- 

vertible monomial changes or positive scalings of variables 
(thanks to Proposition 12. 4p . Also, note that sign(/B(x)J = 

sign (x v Jb{x)) for any x£W + ~ 2 and ugl' _2 . So it clearly 
suffices to show that g(x) := 1 + Xi + ■ • • + xy _ 2 — ix a 
attains only nonnegative values, where now a G R^ -2 and 
ai + - ■ -+aj/_ 2 < 1 (since a,/ lies in the interior of Newt(/s)), 

and 7 = 1' ^ — - (since Condition (b) 

(l-ai aj .,„ 2 )a° 1 ■■■a jl 3 _ 2 



holds with equality) . That g attains only nonnegative values 
is then clearly equivalent to the inequality 

1 +xi H + Xf-2 > jx a 

which is in turn equivalent to 

(1 — ai Oj>_ 2 )(l + xi H h Xj/_ 2 ) 

(1 - ai ay_ a )a; 



> 



n 

i=l 



i'-a 

(1 — CKl Qj'-2) + OilUl H + Qij-'_2%'-2 > f7 u f * 

i=l 

upon substituting Xi = aiUi/(l — a\ — ■ ■ ■ — a.ji-%). The 
last inequality is simply the weighted Arithmetic-Geometric 
Inequality HLP88 so we are done. ■ 



